Coalescing Neutron Stars as Gamma Ray Bursters ? 
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We investigate the dynamics and evolution of coalescing neutron stars. The three- 
dimensional Newtonian equations of hydrodynamics are integrated by the 'Piecewise 
Parabolic Method' (PPM; Colella & Woodward, 1984). The grid is Cartesian and 
equidistant with a resolution of 64 3 or 128 3 , which allows us to calculate the self-gravity 
via fast Fourier transforms. Although the code is purely Newtonian, we do include the 
effects of the emission of gravitational waves on the hydrodynamics in a way originally 
proposed by Blanchet, Damour & Schafer (1990). In particular we treat the gravita- 
tional backreaction in the form found in Shibata et al (1992). These imply solving two 
Q\ additional Poisson equations, which is computationally easy due to the use of FFTs. 

O^j The properties of neutron star matter are described by the equation of state of Lat- 

timer & Swesty (1991), which we handle in tabular form. In addition to the fundamental 
hydrodynamic quantities, density, momentum, and energy, we follow the time evolution 
of the electron density in the stellar gas. Energy loss and changes of the electron abun- 
\q dance due to the emission of neutrinos are taken into account by an elaborate "neutrino 

leakage scheme" , which employs a careful calculation of the lepton number and energy 
t>* source terms of all neutrino types (Janka, 1991; Ruffert, Janka, and Schafer, in prepara- 

^ tion). Neutrinos are produced via thermal processes and lepton captures onto baryons. 

CO Matter is rendered optically thick to neutrinos due to neutrino-nucleon scattering and 

^ absorption of neutrinos onto baryons. 

Q\ We simulate the coalescence of two identical, cold neutron stars with a baryonic mass 

*c3 of r; 1.6 M Q and a radius of r; 15 km. Initially the stars are spherical and orbit around the 

Ph common center of mass with a center-to-center distance of 42 km and an orbital period 

of about 2 ms. The initial distributions of density and electron concentration are given 
^ from a model of a cold neutron star in hydrostatic equilibrium. The temperature is then 

ci chosen in such a way that the thermal energy is set to roughly 3% of the degeneracy level. 

The initial velocity distribution is varied between the three models which we simulate. 
The orbital velocities of the coalescing neutron stars are prescribed according to the 
motions of point masses as computed from the quadrupole formula. Spins are added to 
the neutron stars to take into account rotations around their axes vertical to the orbital 
plane. These spins are different from model to model to describe the cases with spin 
vectors parallel and antiparallel to the vector of the orbital angular momentum, and the 
case with rigid rotation, respectively. 

The orbit decays due to gravitational wave emission, and after one revolution the 
stars are so close that dynamical instability sets in. Within 1 ms they merge into a 
rapidly spinning (P ^ 1 ms), high-density body (p Ri 10 14 g/cm 3 ) with a surrounding 
thick disk of material with densities p r; 10 10 — 10 12 g/cm 3 and rotational velocities 
of 0.3-0.5 c. Because of its high mass of r; 3 M we expect the central object to collapse 
into a black hole essentially instantaneously. 
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The use of a "realistic" EOS in our simulations allows us to determine the temperature 
distribution in the merging stars and means that bulk viscosity effects of the matter are 
taken into account. Initially the neutron stars are cool and at a separation of several 
neutron star radii tidal deformations are insignificant. The neutron stars do not heat up 
noticeably even during the final stages of inspiral, since shear viscosity is small and the 
tidal deformations tend to increase the volume. The parts of the neutron stars closest 
to each other therefore experience tidal stretching and show cooling due to expansion 
rather than heating. Not earlier as when the stars start to merge, shearing effects of 
matter with opposite velocity components and local compression lead to a temperature 
increase. After merging we find peak temperatures of 30-40 MeV in two distinct "spots" 
inside the central, massive object, while the surrounding disk gas (which is only a few 
tenths of a solar mass) remains at temperatures of 4-8 MeV. 

The peak emission of gravitational waves is short but powerful. A peak luminosity 
of a few times 10 55 erg/s is reached for about 1 ms. The amplitudes of the gravitational 
waves reach 3 • 10~ 23 for an assumed distance of 1 Gpc, and the typical frequency is near 
the dynamical value, 2 kHz. 

Combining information about the fluxes of neutrinos and antineutrinos radiated into 
different directions, we are able to evaluate our models for the energy deposition by 
annihilation of vv into e + e~-pairs in the vicinity of the merger. This mechanism was 
suggested as a possible energy source for 7-ray bursts (see e.g., Meszaros & Rees, 1992; 
Mochkovitch et al, 1993). Although the final, dynamical phase of the coalescence happens 
within approximately 3 ms, the individual neutrino luminosities saturate at values of 
several 10 52 erg/s only on the longer timescale of about 8 ms which the outer layers of 
the merged object need to reach their maximum temperatures (partly by shock heating). 
The integral rate of energy deposition by neutrino annihilation gets to a few 10 50 erg/s 
during the phase of strongest neutrino emission. This rate has to be compared with the 
typical energy of 10 50 -10 51 erg observed in 7-ray wavelengths from bursts at cosmological 
distances. In order to yield this energy the merged object would have to radiate neutrinos 
with high luminosities for a period of 1-10 s. This is difficult to conceive in the considered 
scenario, since the merged object with a mass of around 3 M Q will not escape the collapse 
to a black hole on a dynamical timescale of ~ 1 ms. If one takes into account gravitational 
redshift effects and an efficiency of less than 100% for the conversion of e + e~-pair energy 
into 7-rays, our result is off from the desired value by more than a factor of 100. We 
therefore conclude that some other mechanism than ^^-annihilation has to be considered 
if 7-ray bursters are to be identified with coalescing neutron stars. 
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